#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <stdlib.h>
#include <string>
#include <cmath>
#include <time.h>
#include <fstream>
#include <vector>
#include <algorithm> 

using namespace std;



int main() 
{

ifstream infile("peaks_050.txt");
ofstream outfile("scarto.txt");
ofstream outfile2("scarto2.txt");

double t,u, A,dt,dt_min,dt_min_tot,A_min_tot;
int i,j,n,N,v,j_min;


dt_min_tot=0;
A_min_tot=0;
v=0;

N=16000000;

vector<double>t_v,A_v,dt_v;
n=0;
while(infile>>t>>A)
{
t_v.push_back(t);
A_v.push_back(A);
dt_v.push_back(0);
n++;

//cout<<t<<endl;
}

for(int i=0; i< N;i++)
{
u=rand()/(1.0+RAND_MAX);
u=-180.0*u;

for(j=0;j<n;j++)
{
dt_v[j]=sqrt((t_v[j]-u)*(t_v[j]-u));
}
dt_min=1E+06;
for(j=0;j<n;j++)
{
if(dt_v[j]<dt_min){dt_min=dt_v[j], j_min=j;}
}
if(dt_min<1E+05){A_min_tot+=A_v[j_min];dt_min_tot+=dt_min;outfile<<u<<"\t"<<dt_min<<endl;}
if(dt_min<0.5){v++;}

if(i%16==0){outfile2<<u<<"\t"<<dt_min_tot/16.0<<"\t"<<A_min_tot/16.0<<endl;dt_min_tot=0;A_min_tot=0;}
}
cout<<"dt_min_tot:="<<dt_min_tot<<endl;

cout<<"dt_min_tot/N:="<<1.0*dt_min_tot/N<<endl;
cout<<"value:"<<v<<endl;
cout<<"v/N:="<< 1.0*v/N<<endl;




return 0;

}


